Excitation mapping of whispering gallery modes in silica microcavities 
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We report the direct observation of the electromagnetic-field distribution of whispering-gallery modes in silica 
microcavities (spheres and toroids). It is revealed by their excitation efficiency with a tapered fiber coupler 
swept along the meridian. The originality of this method lies in the use of the coupler itself for the near field 
mapping, eliminating the need of additional tools used in previous work. This method is successfully applied 
to microspheres and microtoroids. © 2009 Optical Society of America 
OCIS codes: 140.3945, 060.2430, 180.4243, 130.3990 



Silica optical microcavities formed by surface tension, 
like microspheres [1] and more recently microtoroids [2], 
feature whispering gallery modes (WGM) that have 
drawn a large interest in the last two decades. WGMs are 
characterized by ultra-high quality factors (up to 10 10 
in the infrared) and moderate mode volumes (down to 
about 100 A 3 for a typical diameter of 50 |im) . They have 
led to many applications in various fields like CQED, mi- 
crolasers, telecommunications and biosensing. A WGM 
can be seen as a nearly grazing ray guided around the 
resonator by successive total internal reflections and ful- 
filling a phase matching condition after one roundtrip. 
The "fundamental mode", corresponding to a ray tightly 
bound close to the equator of the cavity, is characterized 
by a single antinode along both the polar and radial di- 
rections, and achieves the smallest mode volume. Most 
applications depending on the mode volume can be opti- 
mized by selectively working with this mode, which needs 
to be unambiguously identified. 

Several papers [3-6] have been devoted to this ques- 
tion. In [3], the near-field of a microsphere is imaged on 
a camera through a coupling prism and a microscope. 
Refs. [4] and [5] are based on direct detection of the 
near-field: a fiber tip is scanned along the sphere surface 
in order to map the evanescent field, for a fixed excita- 
tion frequency. In [6], the near-field of poorly confined 
leaky modes was directly imaged on a camera. The de- 
pendance of the coupling efficiency with respect to the 
mode order that is at the heart of the present paper has 
been used in [7] to filter out high order WGMs of a cm- 
sized microdisk using an auxiliary coupling prism. 

In this paper, we show that near-field mapping does 
not require the use a fiber tip or coupling prism but 
can readily bee obtained using the tapered fiber coupler 
technique [8, 9] and a widely tunable laser. The near- 
field distribution is revealed through the spatial depen- 
dance of the excitation efficiency. We first demonstrate 



this method on microspheres which have a well known 
WGM spectrum and then apply it to a microtoroid, for 
which no simple analytic description exists. 

Beside their polarisation, the WGMs of a microsphere 
are characterized by three integer orders n, I and m. 
The radial order n given by the number of antinodes of 
the radial function is not considered here. The orders 
£ and m are the usual angular momentum quantifica- 
tion numbers, with m G [—•£,•£], corresponding to a field 
distribution approximately proportional to the spherical 
harmonic Y™. The fundamental mode then corresponds 
to the WGM of highest azimuthal number m = £, which 
has a single antinode in the polar direction. It is well 
known [10] that a small ellipticity e breaks the spherical 
symmetry and leads to a frequency shift: 
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This gives a multiplet of nearly equally spaced modes 
for |m| w I with an interval Al> = e v^ m jl m e FSR 
where FSR = c/2irnsa ~ 1000 GHz is the free spectral 




Fig. 1. Sketch of the experimental setup, with definition of 
the x,y,z axes used hereafter. S: Sphere (not to scale), T: 
Taper, PZT: Three-axis piezoelectric stage, G: Goniometer, 
ECLD: Extended cavity laser diode, PD: Photodiode. 
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spectra are displayed as a waterfall in Fig. [3] for a sphere 
of diameter 2a = 56 urn and ellipticity e = 1.2%. The 
successive curves are horizontally offset according to the 
corresponding z value, so that the abscissa corresponds 
to both the coupling efficiency and the z-coordinate. The 
overlapping absorption dips allow a clear visualization of 
the field distribution of the different modes, with q+1 = 
i—\m\ + l antinodes. 
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Fig. 2. Typical taper transmission spectrum for a prolate ^ 30 

microsphere of diameter D ~ 68 |im and ellipticity e ~ 0.4%. O 

Each line corresponds to the q — i — \m\ value given below it. >> 

The bottom curve is the spectrum of a confocal Fabry-Perot § 

providing the frequency scale. a 1 



range of the cavity (a: sphere radius, no,: refraction index 
of silica). It moreover indicates that the resonance fre- 
quency decreases with q = £— \m\ for a prolate spheroid 
(e > 0). Hence the fundamental mode has the largest 
frequency of the n, t multiplet in our experiments with 
prolate microspheres. 

The experimental setup is sketched in Fig. [TJ The mi- 
crosphere is prepared by CO2 laser melting of a tapered 
fiber, allowing to finely control its size and eccentric- 
ity, and to optimize the parallelism of its revolution axis 
with the mother-fiber axis [1]. The taper coupler is pro- 
duced by pulling a single-mode optical fiber heated by a 
microtorch, as described in [11]. A widely tunable mode- 
hop-free extended cavity laser diode (A ~ 775 nm) is 
coupled into the tapered fiber. When scanning the laser 
frequency the WGM resonances are detected by an ab- 
sorption dip on the laser output intensity monitored by a 
photodiode (PD) and a digital oscilloscope. Our setup is 
designed to maximize the excitation of low angular order 
q WGM. The rotation angles around the x and y axes 
and the position along the y- and z-axis should first be 
optimized with the goniometer and PZT shown in Fig. 1. 
The air gap g between the sphere and the coupler is then 
adjusted with the PZT actuator. The gap is set such that 
the resonances are in the undercoupled regime, where the 
broadening of the resonance linewidth is negligible. 

Fig. [2] shows an optimized taper transmission spec- 
trum exhibiting a WGM multiplet. The equally spaced 
lines correspond to the splitting described by Eq. [T] The 
missing resonance on the high frequency side at the loca- 
tion pointed by the dotted line ascertain the identifica- 
tion of the q — WGM, and the subsequent determina- 
tion of the other q values in this series. The "excitation" 
or "coupling efficiency" is defined as the relative depth 
of this transmission dip C = 1 — T res /2o, where Tq is 
the out-of-resonance transmission and T res < Tq is the 
on-resonance transmission. 

We have recorded the WGMs' transmission spectrum 
for different values of the sphere height z. The resulting 
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Fig. 3. Waterfall of transmission spectra obtained for differ- 
ent z-positions: each transmission spectrum is plotted at the 
corresponding z value. 
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Fig. 4. Symbols: coupling efficiency as a function of z- 
position, for q — • • • 5, deduced from Fig. [3] Solid lines: 
global fit of the experimental data. 

For a detailed analysis of these data we use a simplified 
expression of the coupling efficiency valid in the thin 
taper limit, assuming that the overlap integral of the 
WGM and the coupler fields is simply proportional to the 
WGM field. In this model, using spherical coordinates 
(r, 9, 4>) with origin on the sphere center, the coupling 
efficiency writes: 

(o) 

C,(r,0) = K q \ \Yf-\0)\ 2 exp(-2«(r-a)) (2) 

lq 

where K q is a scaling coefficient depending on the ta- 
per diameter and its effective index, 7g ' is the intrin- 
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sic line-width and 7g L ' the observed loaded linewidth, 
k" 1 f» (rig — l)" 1 / 2 A/27T is the evanescent wave char- 
acteristic depth. The ratio of the two linewidths allows 
to take into account overcoupling, that is visible on the 
fundamental mode. 

When substituting the spherical harmonic by its 
Hermite-Gauss asymptotic expression for large I, and 
using a second order expansion of the gap g = r — a = 
go + z 2 /2a, equation J2]) leads to: 



exp 



(3) 



From the experimental results, we have extracted the 
coupling efficiency and linewidth of all the resonances, 
and plotted in Fig. [4] the normalized area of the res- 
onances, given by the product C q {z) x jq°^/^q L \ The 
curves obtained for different q values are offset to evi- 
dence the similarity with Fig. [H 

We then perform a global fit of these curves according 
to Eq. [31 using different amplitudes but the same hori- 
zontal scale for the 6 q values. This fit is plotted as solid 
lines in Fig. [4] and is in very good agreement with the 
experimental points, thus proving that our model is ac- 
curate and that the measured profiles are actually related 
to the field distribution. A departure from the model are 
nevertheless observed: the fitted horizontal z-scale does 
not match the expected one but is larger by a factor of 
about 1.4. This effect does not come from mechanical 
effects or from the PZT calibration, which can easily be 
ruled out. It likely results from the non-negligible taper 
diameter, that is significantly larger than the depth of 
the evanescent wave This idea has been verified in 
a simple numerical simulation of the overlap integral, 
where the finite taper diameter systematically results in 
an effective z smaller than the taper center position used 
in the experiment. 



for their mode structure and resonance positions is not 
available, in particular because the n and I no longer 
exist while m (angular momentum) and q keep their sig- 
nificance. An experimental determination of their field 
distribution is highly desirable. Our new method can 
conveniently be adapted to monitor this field distribu- 
tion. Because of the large curvature of the toroid's minor- 
diameter, a vertical scanning would result in a very large 
gap, and a vanishing signal. Therefore we have replaced 
the linear z scanning by a circular 9 scanning, intended to 
keep an almost constant coupling gap (9 is defined on the 
inset of Fig. [5]). We have used rather thick microtoroids 
(minor diameter d sw 6.5 |ixn), in order to observe several 
modes in a typical frequency range of 50 GHz. The result 
of this experiment is shown in Fig. [5j It is displayed as 
a waterfall, where successive spectra are offset accord- 
ing to 9, from — ir/2 (below the toroid) to ir/2 (above 
the toroid) . This figure shows several modes with clearly 
distinct angular intensity distributions, and in particu- 
lar a fundamental mode at a frequency offset of 19 GHz. 
One can notice that the resonances are broader and more 
pronounced on the left side of the figure because of an 
imperfect centering of the circular motion, leading to a 
smaller gap on this side than on the opposite one, but 
this does not prevent the WGMs' characterization. 

To conclude, we have demonstrated a new and robust 
method to characterize the angular structure of WGMs 
of a microcavity. Based on the taper used for excitation, 
it eliminates steric problem arising for other methods. 
It allows to accurately position the coupling device at 
the equator location, thus optimizing the coupling to the 
most confined mode and canceling odd-modes coupling 
efficiency. 

This work has been supported by C'Nano-Ile de France 
and French embassy in Beijing. G. L. acknowledges sup- 
port from the Chinese Scholarship Council. 
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Fig. 5. Circular taper mapping of a microtoroid. The spectra 
of the transmitted power are offset to the corresponding 6 
position. Inset : definition of coordinate. 

The WGMs of toroidal microcavities have the advan- 
tage to achieve smaller mode volumes than microspheres 
of the same outer diameter. A general analytical theory 
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